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Wc propose a new method for the study of the chiral properties of the ground state in Quantum Field Theories 
(QFT's) which is based on the computation of the probability distribution function (p.d.f.) of the chiral condensate 
in the chiral limit. We show how despite of the fact that Grassmann variables cannot be simulated on a computer, 
an analysis of spontaneous symmetry breaking without a symmetry breaking external field, which is standard in 
the case of spin systems, can also be done in QFT's with fermion degrees of freedom. 



Chiral-symmctry breaking (xSB) is an impor- 
tant feature of many systems with fermions. In 
QCD at finite temperature, there is a chiral phase 
transition from a confining phase, where chiral 
symmetry is spontaneously broken to the de- 
confining phase, where chiral symmetry is re- 
stored. In non-compact QED, a chiral transition 
separates the broken phase and Coulomb phase, 
raising the possibility of the existence of the con- 
tinuum limit. 

Although these issues have been extensively in- 
vestigated for a long time, there are still many 
open questions. The main question we would ad- 
dress here is how to get quantitative information 
on spontaneous chiral-symmetry breaking from 
numerical simulations. 

In the chiral limit, if we use the chiral sym- 
metric action and boundary conditions satisfying 
the symmetry requirement, the value of any or- 
der parameter is always zero, and we don't know 
in which phase the system is. The usual way is 
to add to the action a symmetry breaking term 
so that the fermion condensate is non-zero. An 
extrapolation to the chiral limit has to be done 
after the thermodynamical limit. Unfortunately, 
such an extrapolation in numerical analysis is an 
extrapolation from data obtained on a finite lat- 
tice at small set of fermion masses to the mass- 
less limit, in which the operation is ambiguous 
and therefore the result might be arbitrary. Can 



wc avoid such an ambiguous extrapolation and 
obtain information directly from the symmetric 
action? 

It is well known for bosons or a spin system [1] , 
one can work directly with a symmetry action. 
From the data of a simulation, one can construct 
the probability distribution function (p.d.f.) of 
the order parameter from which the value of the 
order parameter can be estimated by finite size 
scahng [2]. 

For fermions, this method is not straightfor- 
ward. The Grassmann variables can not directly 
be simulated on the computer and Mathews- 
Salam formula has to be used to integrate out 
the fermionic degrees of freedom. Due to such 
an integration, in the chiral limit, the expecta- 
tion value for any order parameter vanishes iden- 
tically. Then it is not possible to construct the 
p.d.f of the order parameter as a histogram. But 
this doesn't mean we can not obtain information 
on spontaneous xSB (SxSB) from the p.d.f. 

First of all, let's introduce some basic concepts 
in our approach. When a symmetry is sponta- 
neously broken, a set of degenerate pure vacua ap- 
pears, forming a Gibbs state, i.e., a macroscopic 
state. The expectation value for any operator O 
in the Gibbs state can be written as 

<0>=^w„<0>„, (1) 
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where a labels the pure vacuum, and Wa is the 
probability of the pure vacuum a in the Gibbs 
state. Let us choose the chiral condensate tptl) as 
the order parameter, then the expectation value 
of the order parameter in the a vacuum is 

^" = ^ 5] < >« . (2) 

All pure vacua, and only pure vacua, satisfy the 
cluster properties, which implies that the inten- 
sive quantities like the magnetization or chiral 
condensate do not fluctuate in the pure vacua. 

With these basic concepts, we can define the 
"macroscopic" p.d.f. of the order parameter c by 

P{c) = Y,w^5{c-c^). (3) 

a 

The function P{c) tells us what is the probabil- 
ity that choosing RANDOMLY a vacuum state, 
we get the value c for the chiral order parame- 
ter. If the vacuum state is invariant under chiral 
transformations, i.e., if it is unique as concern- 
ing the chiral symmetry, P{c) will be a single S 
function S{c). If there is SxSB, P(c) will be a 
complicate distribution. 

Let's go further to look at the form of P{c). 
Since we are interested in the analysis of the chi- 
ral symmetry on the lattice, we will use the stag- 
gered fermion regularization. In such a case, there 
is a continuous U(l) chiral symmetry in the ac- 
tion. Then if this continuous residual symmetry 
is spontaneously broken, we will get a continu- 
ous set of equilibrium states labeled by an angle 
a G [0, 27r]. The v.e.v. of the chiral condensate 
Ca at each a vacuum can be parameterized as 
Co cos (a), Co being the value corresponding to the 
a — vacuum selected when switching-on an ex- 
ternal "magnetic" field and taking the massless 
limit afterwards. In another word, Cq is the chiral 
order parameter in the conventional sense. No- 
tice that no matter in what phase the system is, 
the integration of the fermionic degrees of free- 
dom (equivalently the integration over the pure 
vacua) always leads to < c >= 0. The function 
P{c) can be computed as 



P(c) 
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Figure 1. Standard form of P{c) of (4) in the 
broken phase. 
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P{c) = — / da5{c — cqcos^o)). (4) 
Jo 

In the symmetric phase co = 0, it is obvious 
that P(c) = 5{c). In the broken phase, P{c) = 
l/{7r(cg - c2)i/2) for -Co < c < Co, and P(c) = 
otherwise (see Fig. 1). Its Fourier transformed 

= ^ r d0e">'<^'°'' = Mqco). (5) 
2"" Jo 

Let us emphasis the power of p.d.f.: although 
< c >= in all cases, from the shape of p.d.f., we 
can still determine in which phase the system is 
and what value of the order parameter co would 
be. 

On a finite lattice, let us define p.d.f. as 
Py(c)=<5(l^V;(z)V(x)-c)> (6) 

where the expectation value in (6) is computed 
in the Gibbs state and the integration measure is 
that associated to the partition function 

Z = [ ##d[/e-^«(^)+*^'^. (7) 
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Sg in (7) is the pure gauge action and A the 
fermionic matrix. 

Using the cluster properties of the pure vacua, 
it can be shown [3] that the definition (6) is iden- 
tical to (3) in the thcrmodynamical limit, i.e.. 



limv^ooPvic) = Pic). 



(8) 



Expression (6) is not suitable for numerical 
computation. However its Fourier transformed 
P{q) as will be shown, can be numerically com- 
puted. The fermion matrix A can be decomposed 
as A = m + iA where m is the fermion mass 
and A a hermitian matrix which depends on the 
gauge field configuration. The eigenvalues of A 
are real and symmetric. Taking into account all 
these properties of A, in the chiral limit the fol- 
lowing expression for Pv{q) can be derived 
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Figure 2. Chiral order parameter in the compact 
Schwinger model in 32^ and 64^ lattices at = 
6.344. The solid line corresponds to its continuum 
analytical value. 



Mo) =<X{{1 - ^) > (9) 

3 ^ 

where the product in (9) runs over all positive 
eigenvalues Aj and the mean values arc computed 
with the probability distribution function of the 
effective gauge theory obtained after integrating 
out the fermion fields. The function Pv(q) can 
be computed numerically and then, by inverse 
Fourier transform we get Pv{c)- 

Finite size scaling of (9) allow us to get quan- 
titative information on S^SB. If the symmetry 
is not broken, in the thcrmodynamical limit, the 
right hand side approaches 1. Otherwise, accord- 
ing to (5) and (9), 

limv^^ < - j^2y2 ) >= Joiqco). (10) 

This equation also holds for a given background 
gauge field U, from which we obtain for this given 
configuration the relation between the lattice chi- 
ral condensate Co{U) and ith zeroes of Jq 

co(i,t/) = ^. (11) 

According to finite size scaling, if the small eigen- 
values go to zero slower than 1/V, then Co{U) = 



0, i.e., chiral symmetry is restored for this config- 
uration. If they approach zero as 1/V then in the 
thcrmodynamical limit, then we expect S^SB. Of 
course, on finite volume, it is not surprising that 
a-i/iVXi) depends on i. However in broken phase, 
a plateau for various i should appear. The extent 
of the plateau increases with the lattice size. Av- 
eraging Co{i,U) over different configurations, we 
get the value of the chiral condensate. 

These properties are nicely demonstrated in 
Fig. 2 where we plot the lattice chiral conden- 
sate Co{i) against i for the one-fiavor compact 
Schwinger model at = 6.344. The contin- 

uous line in this figure stands for the continuum 
analytical result (< ■0^ > /e « 0.16) times the 
bare coupHng g. Figure 2 tells us that the ex- 
tent of plateau indeed increases with the lattice 
size and approaches the exact value in the contin- 
uum. Similar behaviors have also been observed 
for other g values in the weak coupling region. Al- 
though the Schwinger model is not physically re- 
alistic, it has been widely accepted as a very good 
laboratory to check new proposals. This is well 
known for two reasons: the one-fiavor massless 
Schwinger model in the continuum limit is exactly 
solvable and it shares many interesting properties 
with other more relevant physical models. There- 
fore, the results for the Schwinger model strongly 
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Figure 3. Pv{c) in 4 flavor non-compact QED at Figure 4. Pv{c) in 4 flavor non-compact QED at 

P = 0.200 on 10^ lattice. /3 = 0.237 on 10^ lattice. 



support the reliability of our method. 

Contrary to the asymptotic free theories, the 
continuum limit of non-compact QED is defined 
through the cliiral phase transition at some finite 
(3. The issue of triviality (or non-triviality) of 
the continuum limit is still under debate (for an 
overview or references see [4,5]). The main diffi- 
culty in establishing the nature of the continuum 
limit is in determination of the critical exponents, 
since it is very sensitive to the critical coupling /9c . 
Recently, the measurements of the susceptibilities 
[6] suggest f3c « 0.202 for 4 fiavors. It is very im- 
portant to have more independent checks on it. 
For this purpose, we extend the p.d.f. analysis 
to this model. Only the results for 4 fiavors on 
10" lattice will be presented here. Figure 3 shows 
Pv{c) at 13 = 0.200. A peak away from the origin 
develops (also for Pv{—c) since it is a symmetric 
function). According to the above arguments, at 
this j3, the system is in the broken phase. Pv{c) 
at /3 = 0.237 is shown in Fig. 4. From the distri- 
bution, we expect that the system is in the sym- 
metric phase at this /3. Therefore figures 3 and 4 
give a consistent picture with [6]. Of course, ac- 
cording to (8), computation on different lattices 
seems necessary. 

In conclusion, we have developed a new quanti- 
tative approach to spontaneous chiral-symmetry 
breaking. The main advantage of this method, 
when compared with standard simulations, is 



that we can work directly in the chiral limit 
and therefore no mass extrapolations are needed. 
From the probability distribution function of the 
order parameter and finite size scaling, we can ob- 
tain not only the information on SxSB, but also 
the value of the chiral condensate. We have tested 
this method in the Schwinger model and apphed 
it to QED. These results should stimulate people 
working in this field to apply this formalism to 
more interesting physical systems, like QCD. 
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